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Abstract 



Contrary to a recent claim, the inferred primordial abundances of the light 
elements are quite consistent with the expectations from standard big bang 
nucleosynthesis when attention is restricted to direct observations rather than 
results from chemical evolution models. The number of light neutrino (or 
equivalent particle) species (N v ) can be as high as 4.53 if the nucleon-to- 
photon ratio (77) is at its lower limit of 1.65 x 1CP 10 , as constrained by the 
upper bound on the deuterium abundance in high redshift quasar absorption 
systems. Alternatively, with N u = 3, 77 can be as high as 8.90 x 10 -10 if the 
deuterium abundance is bounded from below by its interstellar value. 
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In a recent Letter, Hata et al. |J have made the startling claim that the number of 
light neutrino species deduced from considerations of big bang nucleosynthesis (BBN) is 
N v = 2.1 ± 0.3 (ler), i.e. inconsistent with the standard model (iV„ = 3) at the 98.6% 
C.L. Their analysis is based on 3 key inputs. First, they adopt a primordial 4 He mass 
fraction, Y p = 0.232 ± 0.003 (stat) ± 0.005 (syst), estimated from selected observations 
@ of low metallicity HII regions in blue compact galaxies (BCG). Second, they adopt for 
the primordial abundances of D and 3 He, %j2 P = D/H = 3.5l^g x 10~ 5 (95% C.L.) and 
y 3p = 3 He/H = 1.2 ± 0.3 x 10~ 5 (95%C.L.), using a "generic" chemical evolution model |4| 
normalized to solar system abundances and convolved with BBN predictions. Finally, 
they estimate the primordial abundance of 7 Li to be y-j v = 7 Li/H = l^j^x 10~ 10 (95% C.L.), 
and state that this is consistent with observational data on Pop II stars, taking into account 
possible post big bang production and stellar depletion. They now compare all four ele- 
ments simultaneously against the theoretical predictions, with the theoretical uncertainties 
determined by Monte Carlo methods ||, and obtain the likelihood function for N v . This 
yields the best fit quoted earlier, corresponding to an upper bound of N u < 2.6 (95% C.L.). 
Hata et al. suggest that this "crisis" can be resolved if the 4 He mass fraction has been un- 
derestimated by 0.014 ± 0.004 (lcr), or if the constraint on the D abundance is relaxed by 
assuming that the 3 He survival factor g 3 is < 0.10 at 95% C.L. rather than its adopted value 
of 0.25 in their chemical evolution model. 

In an accompanying Letter, Copi et al. 0] perform a Bayesian analysis, adopting the same 
4 He abundance but allowing for a larger statistical error ay, and a possible offset in the 
central value, Y p = 0.232+AF. Departing from previous assumptions ||, they consider three 
different models of galactic chemical evolution which have varying degrees of 3 He destruction 
(with the average survival fraction being about 15% in Model 1, the extreme case), and 
two possible values for the primordial 7 Li abundance estimated from data on Pop II stars, 
7 Li/H = 1.5 ± 0.3 x 10~ 10 (no depletion), and 7 Li/H = 3.0 ± 0.6 x 10~ 10 (depletion by a 
factor of 2). They too conclude that Y p has been systematically underestimated by about 
0.01. Adopting AY = 0.01, allowing for maximal 3 He destruction (Model 1), and assuming 
zero prior for N v < 3, then yields the 95% C.L. bound N„ < 3.5 for the lower 7 Li abundance; 
the bound decreases to 3.4 for the higher 7 Li abundance. Copi et al. |7|] suggest that the 
situation can be clarified by measuring D in quasar absorption systems (QAS) and testing 
whether 7 Li has indeed been depleted in Pop II stars. 

In fact it has been recognized for some time |T(J that the systematic error in Y p may be 



larger than in the above estimate 0]. Recently it has been emphasized |TTJ that previous 



analyses have adopted an old set of helium emissivities |j| while the recalculated new 



set |H| allow much better fits to detailed line ratios. This question has been examined in a 



study |14| of 27 HII regions in 23 BCG identified in the Byurakan surveys, where, as before, 



the regression of the helium abundance against 'metals' (Z) such as oxygen and nitrogen is 
extrapolated down to zero metallicity to extract its primordial value. Whereas use of the 
old emissivities yields Y p = 0.229 ± 0.004 in agreement with previous results, use of the new 
emissivities (together with new correction factors [P| for the collisional enhancement of He I 



emission lines) raises the value to 14 



Y p = 0.243 ± 0.003. (1) 
Both the dispersion of the data points in the regression plots and their slope (dY/dZ « 1.7± 



2 



0.9) are now smaller than found before 0, in better agreement with theoretical expectations. 
The detailed analysis [14] suggests a residual systematic uncertainty in Y p of ±0.001, so we 



adopt the 95% C.L. range 0.236 - 0.250. 

Secondly, D has already been detected through its Lyman-a absorption lines in the 
spectrum of the quasar Q0014+813, due to a foreground QAS at redshift z ~ 3.32 ]TB|. The 
derived abundances lie in the range 

D/H \qas(i) ~ (1-4 — 2.5) x 10~ 4 . (2) 

Further observations have resolved D lines at z — 3.320482 and z = 3.320790, thus elimi- 
nating the possibility of confusion with an 'interloper' hydrogen cloud ||17|| . The measured 



abundances in the two clouds are, respectively, D/H = 1 o~ 3 - 73±0 - 12 and io _3 - 72±0 09 (where the 
errors are not gaussian); an independent lower limit of D/H > 1.3 x 10~ 4 is set on their sum 
from the Lyman limit opacity. Recently, there has been a detection of D/H = l^lgj x 10~ 4 



in another QAS at z = 2.797957 towards the same quasar [|T^]. The errors are higher because 
the D feature is saturated; nevertheless a 95% C.L. lower limit of D/H > 0.7 x 10~ 4 is set. 
There have been other, less definitive, observations of QAS consistent with this abundance, 
e.g. D/H « io-3- 95±0 - 54 at z = 2.89040 towards GC0636+68 @, D/H < 1.5 x 10" 4 at 
z = 4.672 towards BR 1202-0725 |21 and D/H < l(r 3 - 9±0 - 4 at z = 3.08 towards Q0420- 



388 |20|. However, very recently, other observers have found much lower values in QAS at 



z = 3.572 towards Q1937-1009 |22j and at z = 2.504 towards Q1009+2956 f23|]; their average 
abundance is 

D/H \qas(2) = 2.4 ± 0.3 (stat) ± 0.3 (syst) x 10" 5 . (3) 

Unlike the cloud in which the abundance (0) was measured, these QAS also exhibit absorp- 
tion due to carbon and silicon, whose synthesis in stars would have been accompanied by 
destruction of D. It is argued |22| that this must have been negligible since the metallicity 
is very low. Although this is true averaged over the cloud, large fluctuations in the observed 
D abundance are possible since the mass of absorbing gas covering the QSO image is only 
~ 10~ 6 M Q ; thus D may well have been significantly depleted in it by a star which was not 



massive enough to eject 'metals' fl8|| . Keeping in mind that D is always destroyed by stellar 
processing, the high D measurement (Q) in the chemically unevolved cloud should be taken 
as a conservative upper limit on its primordial abundance. (The consequences of assuming 
that this measurement provides the true primordial value have been investigated elsewhere 

More data is clearly needed to establish that there is indeed a "ceiling" to the D abun- 
dance. Nevertheless, these observations already call into question the chemical evolution 
model |4| employed by Hata et al. This model assumes that primordial D is burnt in stars 
to 3 He, a fraction g% > 0.25 of which survives stellar processing when averaged over all stars 
f25| . However, the recent measurement of 3 He/ 4 He = 2.2lQ;g(stat) ± 0.2(syst) x 10~ 4 in 
the local interstellar gas |26] is close to its value of 1.5 ± 0.3 x 10~ 4 in the pre-solar nebula 



||, demonstrating that the 3 He abundance has not increased significantly in the 4.6 x 10 9 
yr since the formation of the solar system. Indeed, g% must be less than 0.1 if an initial 
D abundance as high as in Eq. (^) is to be reduced to its present value in the interstellar 
medium (ISM) [0 
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D/H \ism ~ 1-5 ± 0.2 x 10~ 5 , (4) 

without producing 3 He in excess of its observed value, 3 He/H pa (1 — 4) x 10~ 5 , in galactic 
HII regions ||28|| . This would be so if there is net destruction of 3 He in (1 — 2) M & stars 
through the same mixing process which appears to be needed to explain other observations, 
e.g. the 12 C/ 13 C ratio |29|; a plausible mechanism for this has been suggested recently |30 . 
Note that the ISM abundance of D sets a lower limit to its primordial value y2 P , since there 
is no known astrophysical source of D JJT . 



Finally, accurate 7 Li abundances have been determined for 80 hot, metal-poor Pop II 
stars, of which 3 have no detectable lithium []3"2]] . Ignoring these reveals a trend of increasing 
7 Li abundance with both increasing temperature and increasing metallicity implying that 
about 35% of the 7 Li was produced by galactic cosmic ray spallation processes. The average 



value in the hottest, most metal-poor stars is [32 



7 Li/H \ Pop n = (1.05 - 2.63) x 10~ 10 (95% C.L.). (5) 

However it may be premature to identify this with the primordial abundance until the 
absence of 7 Li in some of these stars is understood. The primordial abundance may instead 
correspond to the much higher value observed in Pop I stars which has been depleted down 
to (and in some cases, below) the Pop II 'plateau'. Stellar evolution modelling |£J then 
indicates that the initial abundance could have been as high as 

7 Li/H \ Pop ! = (9.54 - 15.1) x 10~ 10 (95% C.L.). (6) 



A recent study ||34|| , which includes new data on 7 halo dwarfs, fails to find evidence of 



significant depletion through diffusion, although other mechanisms are not excluded. For 
example, stellar winds can deplete a primordial abundance of 7 Li/H = io _9 - 5±01 down to 
the Pop II value (||) in a manner consistent with observations |3T| . 

In Fig. |l] we show that the standard model (N u = 3) can be consistent with these 
observations over a wide range of the nucleon-to-photon ratio, 77 = n^jn^. It will be close 
to its minimum allowed value if y 2p is actually given by Eq. (0) and y 7p by Eq. (||), while it 
will be close to its maximum allowed value if y 2p is given by Eq. (f|) and y 7p by Eq. (||). Of 
course a value in between is also possible, depending on the true primordial abundances of 
D and 7 Li. Since only the 4 He abundance (Eq. [I]) can be reasonably taken to be primordial, 
we do not attempt a likelihood analysis. Instead we determine the upper bounds on the 
parameters N v and r] corresponding to the two extreme possibilities above, taking into 
account uncertainties in the nuclear cross-sections and the neutron lifetime by Monte Carlo 
methods |6|,|3~6|1. 



First we evaluate how many species of light neutrinos are allowed subject to the upper 
limits Y p < 0.25 (Eq. |), y 2p < 2.5 x 10" 4 (Eq. |) and y 7p < 2.6 x lO" 10 (Eq. |). Fig. | plots 
the number (N) of computation runs (out of 1000) which satisfy the joint observational 
constraints, for different values of N u . It is seen that (up to \/N statistical fluctuations) 
the "95% C.L." bound is N v < 4.53. This bound varies with the adopted upper limit to the 
helium abundance as 

N max = 3 75 + 7g (ym«x _ Q.240). (7) 
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For N„ = 3, Y p should exceed 0.230; to have Y p exceed 0.236 (Eq. |l|), we require rj > 
2.02 x 10~ 10 . 

Secondly, we calculate the maximum value of rj permitted by the data by requiring that 
50 runs out of 1000 (up to ^/N statistical fluctuations) satisfy the constraints y2 P > 1.1 x 10~ 5 
(Eq. § and y lp < 2.6 x 10~ 10 (Eq. §. A good fit for Y p < 0.247 is 

rf™* = [3.19 + 375.7 (Y p max - 0.240)] x lO" 10 ; (8) 

for higher Y p , the yr p constraint does not permit rj to exceed 5.7 x 10~ 10 , similar to the result 
found earlier |36]. If we choose instead to use the more conservative limit y^ p < 1.5 x 10 -9 



(Eq. P), the bound is further relaxed to 



^max 



[3.28 + 216.4 (Y p max - 0.240) + 34521 (F p max - 0.240) 2 ] x 10~ 10 , (9) 



for Y p < 0.252 and saturates at 1.06 x 10~ 9 for higher values, essentially due to the y2 P 
constraint (see Fig. |3|). Thus for Y p max = 0.25, we find rj max = 8.9 x 10~ 10 , which corresponds 
to a nucleon density in ratio to the critical density of = 0.033/i~ 2 , where h is the 
Hubble parameter in units of 100 kms -1 Mpc -1 . For comparison, Hata et al. |TJ find rj = 
4.4tg;g x 10~ 10 (la) and Copi et al. [0] quote the concordance range rj ~ (2 — 6.5) x 10~ 10 . 

In conclusion, the "crisis" is not with big bang nucleosynthesis but with the naive model 
of galactic chemical evolution and restricted set of observational data considered by Hata 
et al. [l|]. Their claim that the standard model is excluded follows from the unreliable 
"95% C.L." limits they adopt on the input elemental abundances. Copi et al. find, using 
statistical arguments, that assumptions concerning primordial abundances in previous work, 
including their own ||, must be relaxed for consistency with the standard model. In fact a 
variety of observational data had already pointed to the need for this, as reviewed elsewhere 

Our results have important implications for physics beyond the standard model ||37||. For 



example, the relaxation of the bound on the number of neutrino species due to the new D 
and 4 He observations now permits the existence of a gauge singlet neutrino even if it has 
large mixing with doublet neutrinos, as in suggested solutions to the solar and atmospheric 



neutrino problems [SS|. On the other hand, if the primordial D abundance is actually close 
to its present interstellar value and the primordial 7 Li abundance is that inferred from Pop I 
stars, then the nucleon density may be consistent with the observed large nucleon fractions 
in clusters of galaxies, even for a critical density universe 
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FIG. 1. Predicted light element abundances versus 77, with 95% C.L. limits determined by 
Monte Carlo. The rectangles indicate the various observational determinations. Only the 4 He 
abundance is established to be primordial in origin. 
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FIG. 2. Number of Monte Carlo runs (out of 1000) which simultaneously satisfy the constraints 
Y p < 0.25, y2 P < 2.5 x 10~ 4 and y-j v < 2.6 x 10 -10 , as a function of -q (in units of 10~ 10 ), for various 
values of N v . 
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FIG. 3. The upper limit to r\ (in units of 10 ) implied by the constraints y2 P > 1.1 x 10 
and y-jp < 1.5 x 10~ 9 , as a function of the maximum value of Y p . 
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